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Magnetic and electronic properties of the recently-discovered material Sr8CaRe3Cu4024 were 
investigated by means of a quantum Monte Carlo simulation, the Green function method and the 
LSDA+[/ (local spin-density approximation plus the Hubbard- [/ term) method. The LSDA+[/ 
calculation shows that the ground state is an insulator with magnetic moment M=1.01/iB/f-u., 
which is consistent with experimental results. The magnetic sites were specified and an effective 
model for the magnetic properties of this compound derived. The resultant effective model is a three- 
dimensional Heisenberg model with spin-alternation. Finite-temperature properties of this effective 
model are investigated by the quantum Monte Carlo method (continuous-time loop algorithm) 
and the Green function method. The numerical results are consistent with experimental results, 
indicating that the model is suitable for this material. Using the analysis of the effective model, 
some predictions for the material are made. 



I. INTRODUCTION 

Strong electronic correlations produce a rich variety 
of phenomena. High-Tc superconductivity 1] is one of 
such a phenomenon. It is widely believed that the ori- 
gin of high-Tc superconductivity is related to the na- 
ture of undoped systems 2], which are Mott insulators 
with an antiferromagnetic order. Recently-discovered 
perovskite cuprate Sr8CaRe3Cu4O24 is a magnetic in- 
sulator which has spontaneous magnetization at room 
temperature and in some respects resembles parent ma- 
terials of high-Tc superconductors. It is expected that 
the ferromagnetism of Sr8CaRe3Cu4024 might be caused 
by strong electronic correlations due to a similar mech- 
anism to the antiferromagnetism of parent materials of 
high-Tc superconductors. Furthermore, among ferromag- 
netic cuprates this material has an unusually high mag- 
netic transition temperature (Tc) up to 440K, although 
those of other ferromagnetic cuprates are at most 30K 4J . 
Thus, the material is interesting for its high Tc induced 
by strong electronic correlations. The purpose of this pa- 
per is to clarify the origin of magnetism and to predict 
the possible magnetic properties of this material. 



II. LOCAL SPIN-DENSITY APPROXIMATION 

The lattice structure is shown in Fig. Q which 
is determined experimentally 3] . In order to investi- 
gate electronic properties of the compound, the local 
spin-density approximation plus the Hubbard- term 
(LSDA-ht/) method using the package WIEN2K 5] which 
is an iniplementation of the density-functional APW-lo 
methodp, has been used. The cutoff values were set 
at i^i^max=7, the LSDA+/7 parameters U=10eV and 
^LSDA+f/=l-2eV 7J. For the exchange correlation, the 
standard generalized gradient approximation (GGA)|^ 
has been applied. The numerical results show that the 
ground-state is an insulator with a ferrimagnetic order: 
The charge gap is estimated as about 1.68eV. Magnetic 
moments are almost localized at Cu sites, and their di- 




FIG. 1: Unit ceh of Sr8CaRe3Cu4024. 



rections at Cul and Cu2 sites are opposite. The abso- 
lute values of the magnetic moments at Cul and Cu2 are 
1.01/iB and 0.86/iB per formula unit (f.u.), respectively, 
and those of other sites are less than 0.07/iB/f-u.. The 
total magnetic moment is l.Ol/iB/f-u., which is compa- 
rable to the experimental result (0.95/iB/f-u.)l3|- The 
LSDA+[/ calculation shows that the orbital degrees of 
freedom are also ordered. The high Tc would be due to 
the large overlap between orbitals of 02 and Cu origi- 
nated from orbital ordering. The mechanism of the an- 
tiferromagnetic coupling is the super-exchange 9j as in 
the parent materials of high-Tc superconductors. Details 
of the LSDA+[/ calculation are presented in a separate 
paper glQj. 



III. EFFECTIVE MODEL FOR MAGNETISM 

By using LSDA+/7 method, the magnetic sites of this 
compound have been identified. Since the system is an 
insulator, it is natural to expect that the effective model 
for magnetism is a Heisenberg model. It was assumed 
that the spins on nearest-neighbor sites interact with each 
other. Then, the Hamiltonian of the effective model is 
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FIG. 2: Unit cell of the effective model, 
expressed as 

n = jY^Si-s,+._, (1) 

where Si and denote spin operators at Cul (i) 

and Cu2 sites, respectively. Here, p represents 

the unit vectors (p=±x, The unit cell of this 

effective model is illustrated in Fig. [3 

Consider the appropriate combination of spin length 
(5, 5), where S and s denote the spin length at Cul and 
Cu2 sites, respectively. In order to determine these, the 
magnetization of the ground state (Mqs) and that of the 
fully-polarized state (Mpp) for the effective model in a 
unit cell were calculated. In the case of (5, 5) = (1/2, 1/2), 
^GS=1 and Mfp=2, which correspond to 2(1^ /i.u. and 
4/iB/f-u., by setting the ^-factor 2. For (6', s)= (1/2, 
1), (1, 1/2), (1, 1), (Mgs, Mfp) becomes (5/2, 7/2), 
(1/2, 5/2), (2, 4), respectively. Here, Mqs and Mfp 
obtained by the LSDA+[/ method are 1.01/iB/f-u. and 
5/iB/f-u.. Hence, (5*, 5) = (1, 1/2) is the most suitable for 
this compound. 

It should be noted that Mqs and Mpp of the effective 
model are exactly 1/2 and 5/2 (equal to Iji^/i.n. and 
5/iB/f-u.) in any size of systems due to the Marshall- 
Lieb-Mattis theorem llj . The ground-state magnetiza- 
tion of the effective model {Iji^/i.n.) is comparable to 
the experimental results 3]. 

Next, the strength of antiferromagnetic coupling J was 
considered. The energy difference between the ground 
state and the fully-polarized state in a unit cell was cal- 
culated by the exact diagonalization method. The re- 
sult is 8J. The energy difference was calculated for 
Sr8CaRe3Cu4024 also by the LSDA+/7 method. The 
result was 0.036i??/. Hence, by comparing these results, 
J=0.0045i??/(=710K) was obtained. 

IV. GREEN FUNCTION METHOD 

Using the above arguments, the effective model 
for magnetism of Sr8CaRe3Cu4024, which is a spin- 
alternating Heisenberg model in three-dimensions has 



been obtained. In order to investigate finite-temperature 
properties of this compound, the Green function 
method [13 was applied to the effective model where 
Green functions are defined as Fourier components of 
time-dependent correlation functions: 

id 

(2) 

where N is the number of unit cells. Since there are 
four sites in a unit cell, 16(=4x4) Green functions are 
necessary. Here, the double-bracket correlation function 
is defined as 

{{A{t);B))^-9mm,B]). (3) 

Coupled equations were obtained for 16 Green functions 
by using a decoupling approximation IT^ 

{{S^t)s+{t);s-))^{S^){s+{t);s-). (4) 

These equations were solved analytically and the temper- 
ature dependence of magnetizations was obtained using 
these Green functions by following Callen's scheme [l^. 
Figure 01 (a) shows the result of the Green function 
method. The agreement with the experimental result 
is good, indicating that the effective model is suitable for 
describing the magnetic properties of Sr8CaRe3Cu4024. 
From the LSDA+/7 result, J was set as J=710K. 

V. QUANTUM MONTE CARLO SIMULATION 

In order to investigate magnetic properties of this 
compound more quantitatively, a quantum Monte Carlo 
method was applied to the effective model. The algo- 
rithm used was the continuous-time loop algorithm (see 
Refs. in 15]), which is a powerful method for non- 
frustrated spin systems. More than one million updates 
for each simulation were performed. The system size was 
up to 16x16x16 unit cells which correspond to 16,384 
sites. 

Spontaneous magnetization can be obtained by extrap- 
olating S{N) which is defined as S{N, T) = yft{S^Sf)^, 
where Sq and Sf are the z-component of Cul-spin at the 
center of the system and that of the furthest site from the 
center, respectively. Here, N denotes the number of unit 
cells. By extrapolating S{N^ T) in the thermodynamic 
limit, we obtain the temperature dependence of sponta- 
neous magnetization as shown in Fig. 01 (b). The quan- 
tum Monte Carlo result is consistent with experimental 
results. By tuning J to fit the transition temperature, we 
obtain J=695K. The first estimation of J=710K based 
on LSDA+/7 results was not far from the fitted value, 
suggesting that the above mapping is valid. To date, 
the temperature dependence of the spontaneous mag- 
netization is the only quantity that has been measured 
experimentally 3] . However, since the effective model for 
this compound has been derived, further predictions of 
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FIG. 3: Temperature dependence of spontaneous magneti- 
zation obtained by (a) the Green function method and (b) 
the quantum Monte Carlo method (sohd curves). Dots and 
dotted hne denote the experimental results 
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FIG. 4: Temperature dependence of the inverse of uniform 
(a) and staggered susceptibilities (b) above the critical tem- 
perature. 



its magnetic properties through the investigation of this 
effective model can be made. 

The temperature dependence of the inverse of uniform 
and staggered susceptibilities above the critical temper- 
ature is shown in Fig. ^ As is evident from this figure, 
both susceptibilities diverge toward the critical point. In 
general, physical quantities show a power-law behavior 
at the critical point of the second-order phase transition. 
The universality class is characterized by the exponents 
of physical quantities such as length scale ^ oc AT~^, 
magnetization M ex |AT|^ and susceptibility x ^ AT~^, 
where AT = T — T^. In finite-size systems, the magne- 
tization and the susceptibility can be expressed in terms 
of scaling functions as 16j 

Mc^ |AT|^M(L|AT|^), (5) 



the three-dimensional (3D) Heisenberg model: i^=0.7054, 
/3=0.3646 and 7=1.3866 17J. The figure shows that data 
in various sizes and various temperatures near the critical 
point fall into single curves, suggesting that the univer- 
sality class is that of the 3D Heisenberg model. 

The temperature dependence of the specific heat is 
plotted in Fig. El The size dependence of the spe- 
cific heat is small. The maximum value is almost sat- 
urated at about 3.1. The critical temperature was esti- 
mated as 0.622±0.002, assuming that the specific heat is 
a maximum at the critical temperature. Since the spe- 
cific heat has little size-effect, it is expected that the spe- 
cific heat at the critical temperature remains finite in the 
thermodynamic limit. This is consistent with the scal- 
ing property of the 3D Heisenberg model: C cx AT~", 
a=-0.1162<0JJJ. 



X-AT ^x(^AT^), (6) VI. CONCLUSIONS 

where L is the number of unit cells in the x-, y- or z- 

direction. M\AT\-^ and xAT^ were plotted as a func- The magnetic and electronic properties of 

tion of L|AT|^ in Fig. O where we use the exponents of Sr8CaRe3Cu4024 have been investigated by means 
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FIG. 5: Scaling plot for M and x- We set the exponents 
of the three-dimensional (3D) Heisenberg model: z/=0.7054, 
/3=0.3646 and 7=1.3866 17]. The data are obtained in clus- 
ters of L=4^16 in the temperature range of T/ J=0.63^0.7. 




FIG. 6: Temperature dependence of the specific heat. Circles, 
squares and diamonds denote the data of L=12, 14 and 16, 
respectively. Solid line is a guide to the eye. 



of the quantum Monte Carlo method, the Green func- 
tion method and the LSDA+[/ method. Since the 
LSDA+/7 calculation shows that magnetic moments are 
localized at Cu sites, a Heisenberg model was introduced 
as an effective model for magnetic properties. The 
spin lengths at Cul and Cu2 sites as one and one half, 
respectively, were obtained by comparing the result 
given by the effective model in a unit cell and that 
given by the LSDA+[/ method. The finite-temperature 
properties of the effective model were investigated by 
the Green function method and the quantum Monte 
Carlo method. The agreement on magnetization at 
finite temperatures between the numerical results and 
the experimental results was good, suggesting that the 
effective model is suitable for describing the magnetic 
properties of Sr8CaRe3Cu4024. The properties of the 
effective model were further investigated and revealed 
that the universality class of the critical point is that of 
the 3D Heisenberg model. Some of the results given by 
the effective model may be accessible by experiments. 
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